First Direct Double-Beta Decay Q- value Measurement of 82 Se in Support of 
Understanding the Nature of the Neutrino 
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In anticipation of results from current and future double-beta decay studies, we report a measure- 
ment resulting in a 82 Se double-beta decay Q-value of 2997.9(3) keV, an order of magnitude more 
precise than the currently accepted value. We also present preliminary results of a calculation of 
the 82 Se neutrinoless double-beta decay nuclear matrix element that corrects in part for the small 
size of the shell model single-particle space. The results of this work are important for design- 
ing next generation double-beta decay experiments and for the theoretical interpretations of their 
observations. 

PACS numbers: 07.75. +h, 14.60.Pq, 23.40.-s, 23.40.Hc 



The results of recent neutrino oscillation experiments 
indicate that the mass of the neutrino is non-zero [H-[3] • 
The mass hierarchy and the absolute mass scale of the 
neutrino, however, are unknown. Furthermore, the na- 
ture of the neutrino is also unknown; is it a Dirac or 
Majorana particle, i.e. is the neutrino its own antiparti- 
cle? The only known practical method for determining 
the nature of the neutrino is through neutrinoless double- 
beta decay (OV/3/3 decay) measurements Interest in 
double-beta decay (/3/3 decay) has been increasing since 
the laboratory verification of the weak, but allowed, two- 
neutrino double-beta decay (2i//3/3 decay) of 82 Se [5|. In- 
cluding laboratory, geochemical, and radiochemical ex- 
periments, twelve isotopes have been observed to undergo 
2z//3/3 decay: 48 Ca, 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 128 Te, 
130 Te, 136 Xe, 150 Nd, 238 U, and 130 Ba @,[3. With the ex- 
ception of the controversial claim in Ref. [8|, OV/3/3 decay 
has yet to be observed. If experiments succeed in ob- 
serving OV/3/3 decay, we would have evidence that the 
neutrino is a Majorana particle and that conservation of 
total lepton number is violated — a situation forbidden 
by the standard model of particle physics. 

An extensive campaign is currently underway to de- 
velop next-generation experiments to detect OV/3/3 decay 
in a number of candidate isotopes (see Ref. Q for a cur- 
rent review of planned experiments). One such experi- 
ment, SuperNEMO, is expected to provide an increase 
in sensitivity of three orders of magnitude over its pre- 
decessor, NEMO-III, and is projected to reach a half-life 
sensitivity at the 90% confidence level of 1-2 x 10 26 years 
by observing 100-200 kg of 82 Se for five years [jJGlj. 

The defining observable of 0^/3/3 decay is a single peak 
in the electron sum-energy spectrum at the /3/3 decay 



Q-value, Qpp. Hence, it is crucial to have an accurate 
determination of Qpp. The Q-value is also a key param- 
eter required to determine the phase space factor (PSF) 
of the decay. The effective Majorana neutrino mass, to- 
gether with the corresponding PSF and nuclear matrix 
element (NME) for a 0^/3/3 decay candidate provide the 
necessary information to determine the 0^/3/3 decay half- 
life, which is given by: 

Go„(Q^ 5 , Z)\M 0v \ 2 {{m^)/m e ) 2 , (1) 
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where M is the relevant NME, {mpp) is the effective Ma- 
jorana mass of the neutrino, m e is the mass of the elec- 
tron, and Gq v is the PSF for the 0^/3/3 decay, which is 
a function of Qpp 5 and the nuclear charge, Z. Thus, to 
obtain an accurate estimation of the half-life sensitivity 
needed to detect a given (mpp), or conversely, to deter- 
mine (mpp) if the half- life is measured, the NME and 
especially the Q-value need to be known to high preci- 
sion. 

Of all the 0^/3/3 decay candidates currently employed 
in experiments, 82 Se is the only one whose Q-value has 
not been measured directly through high-precision Pen- 
ning trap mass spectrometry (PTMS). PTMS has proven 
itself to be the most precise and accurate method for de- 
termining atomic masses and therefore, Q-values ll|. In 
some cases Qp ^-values determined prior to direct Pen- 
ning trap measurements have been found to be off by 



more than ten keV [12j, [13[. Furthermore, careful cal- 
culations of the NME for 82 Se 0^/3/3 decay [Ii-[l8| differ 
from one another by more than a factor of two. In this 
letter, in anticipation of SuperNEMO and other possible 
experiments with 82 Se, we present the results of a direct 
measurement of Qpp with the Low-Energy Beam and Ion 
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Trap (LEBIT) PTMS facility and results to improve shell 
model calculations of the NME. 

The direct Qpp measurement for 82 Se was carried out 
at the National Superconducting Cyclotron Laboratory 
(NSCL) using LEBIT, a Penning trap mass spectrometer 
facility for high precision mass measurements on isotopes 
produced via projectile fragmentation [lij]. We used the 
plasma ion source of the LEBIT facility to simultane- 
ously produce ions of 82 Se and of the (3(5 decay daughter, 
82 Kr. The source was equipped with a ceramic charge 
holder and was filled with ~ 200 mg of granulated se- 
lenium which was vaporized by the heat from the ion 
source filament. A helium support gas for the source was 
mixed with a small amount of krypton to obtain a similar 
beam current of 82 Kr + to 82 Se + within a factor of three. 
The extracted ion beam was guided through a radio- 
frequency quadrupolc (RFQ) mass filter to suppress the 
strong accompanying helium current before injection into 
a cryogenic RFQ beam cooler and buncher for the cre- 
ation of short low-emittance ion bunches that are sent to 
the LEBIT 9.4 T Penning trap 0. On their path the 
ions were purified further by using a time-of-flight mass 
separation scheme 2l|, allowing only ion species of the 
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same mass-to-charge ratio to be dynamically captured in 
the trap. 

At LEBIT, the cyclotron frequency, v c = qB/2Trm, of 
an ion with mass-to-charge ratio, m/q, in a magnetic 
field, B is measured using the Time-of-Flight (TOF) ion- 
cyclotron resonance detection technique [22j,[23|. First, 
isobaric contaminants are removed from the Penning trap 
by driving them to large radial orbits with a radio fre- 
quency (RF) azimuthal dipolc field. To measure v c , the 
trapped ions are exposed to an azimuthal quadrupolc RF 
field at a frequency vjif near their cyclotron frequency 
with the appropriate RF amplitude and excitation time 
[22I [23( 1 . After ejection from the trap, the ions travel 
through the inhomogeneous section of the magnetic field, 
where the energy of the ions' radial motion is trans- 
ferred to the axial direction and are detected on a 
micro-channel plate (MCP) detector. In resonance, i.e., 
vrf = ^c, the energy pickup of the ions' radial motion 
is maximized and their TOF to the MCP is minimized 
[22| . For a cyclotron frequency determination, this cycle 
of trapping, excitation, ejection, and TOF measurement 
is repeated for different frequencies. This leads to cy- 
clotron resonance curves, as shown in Figure 1, with a 
centroid at v c . 

The measurement process for the determination of 
<5^/3( 82 Se) consisted of alternating cyclotron frequency 
measurements of 82 Kr + and 82 Se + . These measurements 
were performed in a series of four runs. The first run con- 
sisted of measurements using a 500 ms excitation time. 
For increased precision, we utilized a 750 ms excitation 
time for the final three runs. Each TOF resonance was 
the average of 25 to 40 scans over the respective fre- 
quency range with 41 trapping cycles per scan. During 
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FIG. 1. (color online). Time-of-flight cyclotron resonance 
curve for 82 Se 



An excitation time of Trf 



750 ms was 

used to obtain a resolving power of 2 x 10 B . The fit of the 
theoretical line shape to the data is represented by the solid 
line. 



the measurement process the number of ions recorded by 
the MCP was limited to an average of 2 ions per trap- 
ping cycle, corresponding to < 7 ions in the trap at a 
time (assuming 30% detector efficiency). This was done 
to limit the number of contaminant ions produced via 
charge-exchange reactions. Each resonance consisted of 
~ 500 to 3000 detected ions, depending on the number of 
scans per resonance and the beam current from the ion 
source. To determine v c , each resonance was fitted using 
the theoretical line shape described in Ref. [22j . The av- 
erage uncertainty in v c for each resonance was ~ 30 ppb 
(parts per billion). 

For the frequency ratio determination of 82 Kr + to 
82 Se + , drifts in the magnetic field were taken into ac- 
count by linearly interpolating two cyclotron frequency 
measurements of 82 Kr + bracketing a 82 Se + measure- 



ment to obtain ^™*( 82 Kr + ). This interpolated cy- 
clotron frequency was used to obtain the frequency ratio 



R = v l c nt ( S2 Ki+)/v c ( S2 Se + ). The values obtained from 
a total of 110 ratio determinations and their weighted 
average are shown in Figure 2. The difference to the 
reference ratio, i?AME2003 = [m( 82 Kr) — ra e ]/[m( 82 Se) — 
m e ], is plotted using the mass values for 82 Kr and 82 Se 
from the atomic mass evaluation AME2003 pHj . 

In preparation for, and during the measurement pro- 
cess, great care was taken to minimize possible system- 
atic effects. By measuring mass doublets, contributions 
to the measurement uncertainty arising from mass de- 
pendent systematic effects due to frequency shifts, for ex- 
ample caused by field imperfections, are already greatly 
reduced. Nevertheless, prior to the measurements, im- 
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FIG. 2. Difference between the cyclotron frequency ratio of 
82 Kr + and 82 Se + and the ratio obtained from literature mass 
data [25l ] . The solid lines indicate the weighted average and 
the 1<t statistical uncertainty band. 



perfections were carefully minimized following a tuning 
procedure described in Ref. j2(| . 

The effect of a non- linear magnetic field drift is not ac- 
counted for in the data evaluation and is not mitigated by 
using a mass doublet. Therefore the cyclotron frequency 
measurements of 82 Kr + and 82 Se + were alternated with 
a period of no greater than 1 hour. Based on an earlier 
study (27j . this should lead to residual systematic effects 
of the cyclotron frequency ratio no larger than 1 ppb. 
This uncertainty was further minimized by stabilizing the 
pressure in the liquid helium dewar to 10 ppm (parts per 
million) [28| , resulting in an uncertainty well below 1 p pb . 
Simultaneously trapped ions of a different m/q value [29( 
can cause frequency shifts — this effect was minimized 
by verifying that contaminant ions were never present at 
a level exceeding a few percent, low enough to not lead 
to a significant shift at the desired precision. 

The weighted average of the ratios for the individ- 
ual runs and the corresponding statistical uncertainty 
are listed in Table I together with the weighted average 
-Rlebit of all results and the value -RAME2003 obtained 
with mass data from AME2003 [25J . Through evaluation 
of the entire data set, with statistical errors as obtained 
from the fits of theoretical lineshapes to the measured 
cyclotron resonance curves, we determined a Birge ra- 
tio [lo| of 1.27(5). While close to unity, the significant 
deviation indicates the presence of residual systematic ef- 
fects at the 0.8 ppb level not discovered in the individual 
measurements or in the tests for systematic effects per- 
formed. Therefore, to account for these non-statistical 
contributions we multiply the statistical uncertainty af 
the weighted average -Rlebit for all data by the value of 



TABLE I. Average cyclotron frequency ratios Rmn= 
v %nt ^82 Kr + ) ( 82 Se + ) with their statistical errors as obtained 
in four separate runs with N frequency ratio measurements 
performed in each run. Also given is the final weighted aver- 
age -Rlebit with its statistical and final uncertainty and the 
ratio based on the AME2003 atomic mass evaluation 12511. 



Run 


N 


Rrun 


1 


53 


1.000 039 285(5) 


2 


2 


1.000 039 30(2) 


3 


7 


1.000 039 29(1) 


4 


48 


1.000 039 290(4) 


Rlebit 

RaME2003 




1.000 039 290(4) (5) 
1.000 039 26(3) 



the Birge ratio. Both the statistical and total uncertainty 
for -Rlebit are given in the table. 

The /3/3 decay Q-valuc is determined from the mass 
difference between the mother nuclide of mass m m and 
daughter nuclide of mass md through: 



Q/3/3/c 2 



m d 



(R - 1) (m d - m e ), (2) 



where R is the cyclotron frequency ratio between the 
singly charged ions of the daughter and mother nuclides, 
c is the speed of light, and m e accounts for the missing 
electron mass of singly charged ions used in the measure- 
ment. Using our final frequency ratio -Rlebit and the 
AME2003 mass for 82 Kr we obtain, Q pp = 2 997.9(3) 
keV. The new LEBIT Q-value is nearly an order of mag- 
nitude more precise than the previous value of Qpp = 
2 996(2) keV based on mass data from [25j ]. a dramatic 
improvement to one of the ingredients needed for a bet- 
ter determination of the half-life limit for 0i//3j3 decay in 
82 Se+. 

In addition to a precise Qpp-value, an accurate NME 
for the decay is also needed. A number of methods have 
been applied to the NME problem, and in 82 Se one of 
the most prominent methods — the shell model 3 H| 
— gives matrix elements that are less than half the size 
of those produced by the quasiparticle random phase ap- 
proximation (QRPA) the Interacting Boson Model 
[13], and the Generator Coordinate Method [3]. Each 
approach has its strengths and weaknesses; for instance, 
the shell model incorporates complicated correlations but 
only in a small single-particle space (the valence shell), 
while the QRPA is applied in large single-particle spaces 
but includes only simple correlations. Here we report an 
attempt to correct for the shell model deficiency by cal- 
culating an effective 0^/3/3 decay operator that implicitly 
includes effects of single-particle levels that are outside 
the valence shell. 

Our calculation, which will be described in detail 
in a forthcoming publication [3lj ]. uses diagrammatic 
perturbation theory to construct an effective two-body 
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TABLE II. Se Ov/3/3 decay matrix elements calculated in 
QRPA, standard shell model (SM), and with the effective 
Qv/3/3 decay operator discussed in the text. 



QRPA [16] 


SM [15] 


Corrected SM 


5.19 


2.64 


3.56 



the 90% confidence level of 1-2 x 10 26 years, an effective 
neutrino mass as low as 60-85 meV could be detected. 
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Ovpp decay operator that, together with an effective 
Hamiltonian, allows the shell model to get around its 
truncations when the procedure is carried to completion 
32l |33j. This framework has been used extensively to 
determine effective valence-space interactions, but much 
less so to construct other effective operators. The only 
existing work on an effective 0v/3j3 decay operator is an 
exploratory calculation, as it happens, for 82 Se [34[. In 
that work the starting point was a G-matrix interac- 
tion [33| . Using a low- momentum interaction 35 1 , Vj ow k , 
derived from nuclear forces in chiral effective field the- 
ory [36[ , we have carried this calculation significantly fur- 
ther. We now include all diagrams to second order in the 
interaction, state norms, and folds, while expanding the 
set of high- lying single-particle orbits that we treat. 

Our result is shown in Table [Til which we compare 
to values obtained from QRPA and the standard shell 
model. The contributions from outside the valence space 
to the effective operator increase the shell model NME 
by ~ 30% to 3.56, narrowing the gap between it and the 
QRPA. 

This calculation represents an important first step in 
producing a true ab initio NME, suggesting that it will 
be larger than the shell model has indicated. We plan to 
improve our results by including three- nucleon forces 37- 
[39l | and two-body currents , by constructing an effec- 
tive interaction consistent with the effective decay oper- 
ator, and by investigating the size of induced three-body 
terms in the effective decay operator. 

In conclusion, by using Penning trap mass spectrome- 
try we have performed the first direct Q-value measure- 
ment of 82 Se /3/3 decay by measuring the cyclotron fre- 
quency ratio between singly charged ions of 82 Se and the 
PP decay daughter, 82 Kr. Our result, Q p/3 = 2 997.9(3) 
kcV, is nearly an order of magnitude more precise than 
the previous value based on the 2003 atomic mass eval- 
uation [i^]. Following the procedure in Ref. 41 1 and 
using our new Q-value, we calculate the PSF for the 
Ov0P decay mode of 82 Se to be G 0v = 2.848(1) x 10" 14 
yr , where the uncertainty has also been improved by 
nearly an order of magnitude. With the corrected shell 
model NME calculation presented here and the current 
upper limits of {mpp) = 140 - 380 meV from the EXO-200 
experiment [IH we obtain a lower limit range for the 82 Se 
QvP0 decay half-life of 5.0 x 10 24 - 3.7 x 10 25 years. As- 
suming SupcrNEMO achieves its projected sensitivity at 
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